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ABSTRACT:  This  document  describes  spatially  explicit  methodology  to  calibrate  the  Land  use  Evolution  and  impact 
Assessment  Model  (LEAM)  using  historic  land  use  data.  The  LEAM  model  is  used  to  predict  urban  growth  and 
environmental  impacts  associated  with  urban  growth.  The  model  uses  physical  and  social  factors  to  simulate  the 
dynamics  of  urban  development.  The  Procedure  for  Estimating  the  Weights  of  LEAM  (PEWL)  uses  logistic  regression 
techniques  and  historic  land  use  data  to  statistically  calibrate  the  13  sub  factors  used  in  LEAM.  PEWL  consists  of 
computer  programs  to  prepare  samples,  generate  programs  for  regression  analysis,  and  produce  data  files  (estimated 
weights)  for  simulation  using  the  Urban  Sprawl  Submodel  of  LEAM.  The  generated  analysis  programs  from  PEWL  are 
executed  using  SAS,  a  common  statistical  package,  to  complete  logistic  estimation.  PEWL  was  developed  in  FORTRAN 
language  on  Microsoft™  Developer  Studio  (FORTRAN  PowerStation  4.0,  1993-1994)  software.  PEWL  includes  both 
the  source  and  executable  files  of  the  FORTRAN  programs.  For  demonstration  purposes,  the  approach  is  demonstrated 
for  a  multi-county  area  around  Fort  Benning,  GA. 
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1  Introduction 

Background 

The  LEAM  Model 

This  document  describes  and  demonstrates  a  quick  and  inexpensive  approach  to 
calibrate  the  Land  use  Evolution  and  impact  Assessment  Model  (LEAM)  using  his¬ 
toric  land-use  data.  LEAM  is  used  to  predict  urban  growth  and  environmental  im¬ 
pacts  associated  with  urban  growth.  It  uses  physical  and  social  factors  to  simulate 
the  dynamics  of  urban  development. 

The  Procedure  for  Estimating  the  Weights  of  LEAM  (PEWL)  uses  logistic  regression 
techniques  and  historic  land-use  data  to  statistically  calibrate  the  13  subfactors 
used  in  LEAM.  PEWL  consists  of  computer  programs  to  prepare  samples,  generate 
programs  for  regression  analysis,  and  produce  data  files  (estimated  weights)  for 
simulation  using  the  Urban  Sprawl  Submodel  of  LEAM.  The  generated  analysis 
programs  from  PEWL  are  executed  using  SAS,  a  common  statistical  package,  to 
complete  logistic  estimation.  PEWL  was  developed  in  FORTRAN  language  on  Mi¬ 
crosoft®  Developer  Studio  (Fortran  PowerStation  4.0,  1993-1994)  software.  PEWL 
includes  both  the  source  and  executable  files  of  the  FORTRAN  programs.  For  dem¬ 
onstration  purposes,  the  approach  is  demonstrated  for  a  multi-county  area  around 
Fort  Benning,  GA. 

Land-use  Conflicts 

Incompatible  land-use  conflicts  that  affect  military  installation  training  and  testing 
are  developing  as  the  population  of  the  United  States  continues  to  grow  and  com¬ 
munication  and  transportation  technologies  support  more  diffused  settlement  pat¬ 
terns  (Goran  et  al.  2000).  Various  techniques  have  been  developed  and  applied  to 
predict  the  potential  for  incompatible  land  uses  developing  around  military  installa¬ 
tions  (Deal  2002;  Deal  et  al.  2002;  Orland  et  al.  2001).  Techniques  range  from  sim¬ 
ple  questionnaires  answered  by  installation  personnel  to  spatially  explicit  dynamic 
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simulation  modeling.  Perhaps  the  most  common  approach  for  predicting  the  poten- 
tial  for  future  incompatible  land  use  involves  applying  1-,  3-,  and  5-mile  buffers 
around  installations  using  current  and  historic  digital  land-use  maps.  By  counting 
the  urbanized  area  in  these  buffers  over  a  time-series  of  maps,  a  simple  graph  of 
time  vs.  total  amount  of  urbanized  area  can  be  created  (Lozar  et  al.  2003).  These 
simple  trends  can  be  very  useful  and  are  often  sufficient  for  many  analyses. 

The  extrapolation  of  past  trends  can  be  misleading,  however.  Growth  might  con¬ 
tinue  to  accelerate  in  situations  where  land  available  for  development  is  still  plenti¬ 
ful  and  where  the  fringe  of  a  major  urbanized  area  is  entering  the  area  of  interest. 
Growth  might  be  moderate  if  an  urbanized  center  is  already  in  the  buffer  and  grow¬ 
ing  along  the  edge  of  an  installation.  It  might  decelerate  if  little  or  no  land  is  left  to 
develop.  To  help  remove  the  uncertainty,  a  more  careful  analysis  of  the  spatial  rela¬ 
tionships  of  growing  urban  centers,  available  land,  natural  (e.g.,  rivers)  or  man¬ 
made  barriers  (e.g.,  limited  access  highways),  and  zoning  can  be  used  to  predict  po¬ 
tential  incompatible  land-use  situations.  LEAM  was  developed  specifically  for  this 
purpose  (Deal  2002). 

LEAM  is  a  computer-based  model  that  simulates  land-use  change  across  space  and 
time.  It  enables  planners,  policy-makers,  interest  groups,  and  laypersons  to  visual¬ 
ize  and  test  communal  decisions  and  evaluate  potential  consequences  of  those  deci¬ 
sions.  The  LEAM  environment  is  used  to  enhance  our  understanding  of  the  connec¬ 
tion  between  urban,  environmental,  social,  and  economic  systems. 

The  fundamental  LEAM  approach  to  modeling  urban  land-use  transformation  dy¬ 
namics  begins  with  drivers,  those  forces  (typically  human)  that  contribute  to  land- 
use  change.  Each  driver  is  developed  as  a  contextual  submodel  run  simultaneously 
in  each  grid  cell  of  raster-based  geographic  information  system  (GIS)  map(s);  linked 
to  form  the  main  framework  of  the  model  and  produce  landscape  simulation  scenar¬ 
ios.  Submodels  are  completed  and  run  independent  of  the  larger  LEAM  framework 
so  that  variables  can  be  scaled  and  plotted  in  formats  that  help  visualize  and  cali¬ 
brate  submodel  behavior  before  it  becomes  integrated  into  the  larger  model. 

LEAM  model  drivers  represent  the  dynamic  interactions  between  the  urban  system 
and  the  surrounding  landscape.  These  drivers  have  associated  weights  that  are 
used  to  combine  the  drivers  into  a  probability  for  urbanization.  Scenario  maps 
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visually  represent  the  resulting  land-use  changes.  Altering  input  parameters 
(different  policies,  trends,  and  unexpected  events)  change  the  spatial  outcome  of  the 
scenario  being  studied.  This  enables  “what-if’  planning  scenarios  that  can  be 
visually  examined  and  interpreted  for  each  simulation  exercise.  A  critical  com¬ 
ponent  to  using  the  LEAM  model  is  calibration  of  the  weights  for  the  individual 
model  drivers. 


Objective 

The  objective  of  this  project  was  to  create  a  methodology  and  tool  set  to  calibrate 
driver  weights  for  LEAM  that  (1)  uses  readily  available  land-use  data,  (2)  re-quires 
minimal  data  preparation,  (3)  is  easily  implemented,  and  (4)  is  adaptable  for  vari¬ 
ous  implementations  of  LEAM.  The  methodology  must  provide  optimum  driver 
weight  calibrations  based  on  available  calibration  data. 


Approach 

LEAM  is  currently  calibrated  using  subject  matter  expert  opinion.  An  alternate 
calibration  approach  was  identified  as  a  critical  component  of  the  software  develop¬ 
ment  process.  Calibration  methods  used  in  a  prior  modeling  effort  were  identified 
as  having  application  in  the  LEAM  system  (Fang  et  al.  2002).  The  calibration  meth¬ 
ods  were  modified  for  the  LEAM  system  and  compared  against  historic  data  to  ver¬ 
ify  that  the  methods  were  applicable.  Validation  results  are  documented  in  Fang  et 
al.  (2005).  Once  the  approach  was  validated  with  historic  data,  the  methods  were 
automated.  The  resulting  methodology  and  tool  set  is  described  in  detail  in  this  re¬ 
port  and  demonstrated  with  data  for  Fort  Benning,  GA.  All  instructions,  software 
code,  and  scripts  are  included  in  appendices  that  allow  users  to  conduct  the  analysis 
for  their  areas  of  interest  and  modify  as  required  for  future  versions  of  LEAM. 


Scope 

Future  training  and  testing  capacities  on  most  installations  will  be  directly  and  in¬ 
directly  affected  by  regional  land-use  policies  and  investments.  Many  of  the  en¬ 
croachment  factors  affecting  training  are  a  direct  result  of  the  surrounding  regional 
urban  patterns.  These  patterns  are  the  result  of  free-market  property  exchanges 
and  ownership.  Land  attractiveness  to  free  market  participants  (e.g.,  home  owners) 
is  based  on  local  county  and  city  investments  and  policies  including  roads,  utilities, 
and  zoning.  Future  training  and  testing  capacities  are,  in  part,  a  function  of  these 
regional  land-use  investments  and  policies.  LEAM  is  a  tool  set  developed  to  help 
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installations  assess  land-use  policies  and  is  applicable  to  all  installations  that  have 
sufficient  data  to  populate  the  model.  The  effort  described  in  this  technical  report 
provides  model  calibration  tools  to  support  implementation  of  LEAM.  The  efforts 
described  in  this  report  are  applicable  to  all  site-specific  implementations  of  LEAM. 
Specific  LEAM  calibrations  shown  in  this  report  are  for  demonstration  purposes 
only. 


Mode  of  Technology  Transfer 

This  document  and  accompanying  scripts,  code,  and  digital  maps  provides  part  of 
the  technology  transfer  of  the  project  to  military  installations,  their  supporting  or¬ 
ganizations,  and  the  Office  of  Economic  Adjustment. 

This  report  will  be  made  accessible  through  the  World  Wide  Web  (WWW)  at  URL: 
http://www.cecer.army.mil 
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2  Overview  of  Procedures  for  Estimation 
of  the  Weights  of  LEAM 


LEAM  predicts  urban  sprawl  and  its  environmental  impact.  It  is  a  cellular  auto¬ 
mata  model.  One  of  the  core  submodels  of  LEAM  uses  physical  (natural  and  con¬ 
structional)  and  social  factors  to  simulate  the  dynamics  of  urban  development  (Ur¬ 
ban  Sprawl  Submodel).  Previous  studies  developed  methods  to  estimate  the 
weights  (coefficients)  of  the  factors  using  logistic  regression  (see  Appendix  A).  Cur¬ 
rently,  more  than  13  possible  factors  are  incorporated  into  this  submodel.  Due  to 
the  amount  of  available  information  for  a  specific  region,  however,  not  all  the  possi¬ 
ble  factors  may  be  available  for  modeling.  This  situation  requires  that  the  proce¬ 
dure  to  estimate  the  weights  of  the  factors  be  flexible.  PEWL  uses  a  general  logistic 
regression  solution  and  characteristics  of  the  Urban  Sprawl  Submodel.  The  PEWL 
procedure  uses  computer  programs  to  prepare  data  samples,  generate  program 
scripts  for  regression  analysis,  and  produce  data  files  (estimated  weights)  for  simu¬ 
lation  using  the  Urban  Sprawl  Submodel.  The  generated  analysis  programs  from 
PEWL  are  executed  using  SAS,  a  very  common  statistical  package,  to  complete  lo¬ 
gistic  estimation. 

PEWL  was  developed  using  the  FORTRAN  language  in  Microsoft®  Developer  Studio 
Fortran  PowerStation  4.0,  1993-1994)  software.  PEWL  includes  both  the  source 
and  executable  files  of  the  FORTRAN  programs. 
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3  General  Settings 


PEWL  requires  a  general  computer  directory  (folder)  and  a  set  of  subdirectories 
(subfolders).  The  suggested  name  of  the  directory  is  “LEAM_Calb”.  The  executable 
files  (and  their  FORTRAN  source  files)  of  PEWL  should  be  copied  into  this  folder. 
The  guide  file  “Parameters.dat”,  which  contains  the  essential  information  for  the 
estimation  procedure,  should  also  be  stored  in  the  LEAM_Calb  folder.  When  a 
PEWL  user  compiles  the  FORTRAN  source  files,  they  need  to  make  sure  that  the 
generated  executable  files  are  in  this  folder  and  not  the  “debug”  subfolder,  which 
may  be  automatically  generated  by  some  FORTRAN  compilers. 

At  least  three  subfolders  must  be  created  under  the  “LEAM_Calb”  directory.  Sub¬ 
folder  “DataFiles”  is  for  original  input  data  files,  which  include  digital  maps  of  land 
use/cover  and  model  factors  (called  attractors  in  LEAM  documentation).  These  digi- 
tal  maps  must  be  in  ASCII  format.  A  subfolder  named  “Sampled”  is  for  storage  of 
the  data  files  that  contain  the  random  samples  for  estimation.  The  subfolder  “Es- 
timatn”  is  for  SAS  programs  and  their  output  files. 

PEWL  contains  a  guide  file  and  three  executable  (and  FORTRAN  source)  files.  Ta¬ 
ble  1  lists  the  steps  from  setting  up  PEWL  to  generating  the  final  product,  the  esti¬ 
mated  weights  of  LEAM.  The  details  of  each  step  are  described  in  subsequent  chap¬ 
ters. 


ASCII  =  American  Standard  Code  for  Information  Interchange 
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Table  1. 

PEWL  set  up. 

Step 

Task 

Product/Location 

System 

1 

Set  up  PEWL  in  computers 

Copied  files  in  created  folders 

Operation 

2 

Prepare  digital  map  files 

Maps  in  Subfolder  “DataFiles” 

GIS 

3 

Create  the  guide  file 

File  “Parameters.dat”  in  Folder  “LEAM_Calb” 

Notepad 

4 

Draw  random  samples  (run 
“LEAM_Sampling.exe”) 

Data  files  in  Subfolder  “Sampled” 

PEWL 

5 

Generate  SAS  programs  (run 
“  L  E  AM_S  AS_B  u  i  Id  e  r.  exe”) 

SAS  programs  in  Subfolder  “Estimatn” 

PEWL 

6 

Estimate  the  weights  (run  SAS 
programs  generated  at  Step  5) 

Estimates  on  SAS  Output  screen 

SAS 

7 

Create  data  files  of  estimates 

Data  files  of  estimates  in  Subfolder  “Esti¬ 
matn” 

Notepad 

8 

Integrate  the  estimates  (run 
“LEAM_Est_lntegratn.exe”) 

Data  file  “Est_All.dat”  in  Subfolder  “Estimatn” 

PEWL 
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4  Guide  File 


PEWL  uses  a  guide  file  named  “Parameters.dat”.  The  guide  file  contains  all  essen¬ 
tial  parameters  for  drawing  random  samples,  generating  SAS  programs,  and  inte¬ 
grating  the  output  files  of  SAS  programs.  Figure  1  shows  the  format  of  this  file. 

The  first  line  of  the  guide  file  is  the  title  line,  which  is  a  brief  description  of  the  pur¬ 
pose  of  the  guide  file.  The  guide  file  has  seven  sections,  each  of  which  starts  with  an 
indication  line  (information  line).  Information  lines  provide  the  in-formation  for  us¬ 
ers  to  input  correct  parameters  in  the  corresponding  sections.  Both  title  and  infor¬ 
mation  lines  can  be  modified  according  to  the  user’s  preference.  These  lines  cannot 
be  removed,  however,  and  each  information  or  title  line  should  occupy  exactly  one 
row.  No  space  line  (row)  is  allowed  between  any  two  lines  (rows)  of  title/  informa¬ 
tion  and  subsequent  parameters. 

The  guide  file  contains  parameters  for  the  folder  location,  the  dimension  of  all  maps, 
and  the  file  names  of  land-use  and  factor  maps.  The  first  parameter  section  pro¬ 
vides  the  location  of  the  general  folder  “LEAM_Calb.”  The  location  should  be  a 
drive  or  folder  (directory)  named  with  fewer  than  60  characters.  When  inputting 
the  location  of  the  general  folder,  it  is  important  to  key  in  the  “\”  sign  as  the  last 
character,  to  indicate  a  folder  name.  For  example,  when  the  location  of  the  general 
folder  “LEAM_Calb”  is  “D:\test\”  and  the  directory  of  “LEAM_Calb”  is 
“D:\test\LEAM_Calb\”,  the  correct  location  parameter  in  the  guide  file  should  be 
“D:\test”,  instead  of  “D:\test\”. 

Section  2  of  the  guide  file  provides  the  dimension  of  all  digital  maps.  The  first  di¬ 
mension  parameter  is  the  number  of  pixel  rows  of  digital  maps.  The  second  pa¬ 
rameter  is  the  number  of  pixel  columns  of  the  maps.  The  dimensions  of  all  maps 
have  to  be  the  same  before  generating  the  guide  file. 

Section  3  provides  the  file  names  of  land-use/cover  maps.  Two  historical  maps  must 
be  provided  in  chronological  order.  The  first  map  represents  the  oldest  period  of 
time. 

The  next  four  sections  of  the  guide  file  provide  the  information  about  the  factors 
used  in  logistic  modeling.  In  the  example  (see  Figure  1),  a  “factor”  is  called  an  “at¬ 
tractor”  in  the  Information  Lines.  First,  the  number  of  factors  needs  to  be  given 
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(Section  4).  Then  the  file  names  of  the  factor  maps  are  listed  by  the  classes  of  de¬ 
velopment.  Sections  5,  6,  and  7,  respectively,  list  the  file  names  for  modeling  the 
probability  of  possible  land  use  converted  to  Residential,  Commercial-Industrial, 
and  Open  Space.  The  order  of  the  classes  of  development  should  not  be  changed.  In 
each  of  these  sections,  under  the  Information  Line,  the  number  of  file  names  in  the 
list  has  to  be  the  same  as  given  in  Section  3.  As  long  as  a  factor  map  is  used  for 
modeling  the  probability  of  one  kind  of  land-use  conversion,  the  file  name  of  that 
factor  has  to  be  included  in  the  file  list  of  that  kind  of  conversion,  no  matter  if  the 
file  name  of  that  factor  has  or  has  not  been  listed  in  the  sections  for  other  kinds  of 
land-use  conversion.  Therefore,  if  a  factor  has  the  same  effect  on  all  three  kinds  of 
land-use  conversion,  the  same  file  name  of  that  factor  has  to  be  included  in  all  file 
list  sections.  Users  can  determine  the  order  of  the  file  name  lists  according  to  their 
preferences.  The  order  must  be  the  same  for  all  file  name  lists,  however,  since 
PEWL  uses  the  order  of  the  file  name  lists  to  sort/arrange  the  factors.  Thus,  vary¬ 
ing  the  order  of  files  in  different  file  lists  will  cause  modeling  errors. 


The  Parameters  of  Sampling  for  Calibration  of  the  Weights  of  LEAM.  < - File  Title 

Folder  Location:  <(■ -  Information  line  of  Section  1 

C:  <5 - The  folder  is  located  right  below  Driver  C 

Map  Dimension.  <( - Information  line  of  Section  2 

NumberOfRows:  500  < - Number  of  rows 

NumberOfColumns:  500  < - Number  of  columns 

Tire  Names  of  Land  Use  Map  Files:  ^ - Information  line  of  Section  3 

IstMap:  stlme_landcover_1992.asc  < - The  name  of  the  1st  map 

2ndMap:  stlme_landcover_2000.asc  < - The  name  of  the  2nd  map 

The  Number  Of  Attractors:  -  Information  line  of  Section  4 


Tire  List 


The  List 


The  List 


3  - Number  of  attractors  considered,  p=3 

of  Attractor  Files  for  Residential  (RES)  Conversion:  < - 

RoadAttRes.asc  <( - The  file  name  of  the  1st  attractor,  x-j 

SlopeRes.asc  <: - The  file  name  of  the  2nd  attractor,  ^ 

Water AttRes.  as  C  y - The  file  name  of  the  last  attractor,  xp 

of  Attractor  Files  for  Commercial-Industrial  (Cl)  Conversion: 

RoadAttCom.asC  ^ - The  file  name  of  the  1st  attractor,  x  j 

SlopeCom.asc  4 - The  file  name  of  the  2nd  attractor,  X2 

Water  AttCom.  as  C<( - The  file  name  of  the  last  attractor,  Xp 

of  Attractor  Files  for  Open  Space  (OS)  Conversion:  < - 

RoadAttOS.asc  y  The  file  name  of  the  1st  attractor,  x-j 

SlopeOS.asc  4  The  file  name  of  the  2nd  attractor,  X2 

W  aterAttOS.asc  <(  The  file  name  of  the  last  attractor,  xp 


Information  line 
of  Section  5 


^-1 

Information  line 
of  Section  6 

Information  line 
of Section  7 


**********3(:************^£n(j  0f  File************************ 


H 


■  Termination  of  the  guide  file 


Figure  1 .  Format  of  the  guide  file  (example). 
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5  Sampling 


Based  on  parameters  provided  in  the  guide  file,  random  samples  will  be  generated 
by  running  the  program  “LEAM_Sampling.exe”  (or  by  first  compiling 
“LEAM_Sampling.f  ’  to  obtain  the  executable  file).  Before  running  this  program,  all 
digital  maps  named  in  the  guide  file  have  to  be  prepared  and  stored  in  the  subfolder 
(subdirectory)  “DataFiles.” 

Sample  size  is  automatically  determined  by  the  PEWL  program.  When  the  total 
number  of  pixels  in  a  digital  map  is  less  than  350,000,  about  2.5  percent  of  pixels 
will  be  randomly  sampled.  Otherwise,  about  9,500  pixels  will  be  randomly  sampled 
for  estimation.  The  random  samples  are  saved  in  four  data  files  in  the  subfolder 
“Sampled.”  The  names  of  the  data  files  are  “RES_Sampled.dat,”  “CI_Sampled.dat,” 
“OS_Sampled.dat,”  and  “Nch_Sampled.dat,”  which  are  used  respectively  for  esti¬ 
mating  the  weights  (coefficients)  of  probability  models  for  converting  to  Residential, 
Commercial-Industrial,  and  Open  Space  use,  and  to  “No  development.” 

Each  of  the  four  data  files  has  the  same  general  structure.  In  a  data  file,  the  first 
column  is  the  code  of  land-use  conversion.  When  a  sampled  pixel  changed  from 
available  to  a  specific  land  use  (Residential,  Commercial-Industrial,  or  Open  Space) 
as  indicated  by  the  file  name,  the  code  is  0.0000;  otherwise,  the  code  is  1.0000.  This 
code  system  is  based  on  the  definitions  of  SAS  logistic  procedure,  which  uses  0  (zero) 
to  represent  the  positive  response  (happened)  of  an  event.  From  the  second  column 
to  the  last  (p+1)  column,  there  are  totally  p  (the  number  of  factors  used  for  modeling 
as  defined  in  Section  3  of  the  guide  file)  columns.  The  order  of  these  columns  corre¬ 
sponds  to  those  defined  in  the  data  file  list  in  the  guide  file.  For  example,  in  the 
sampled  data  file  “CI_Sampled.dat,”  the  first  column  is  the  code  to  indicate  the  con¬ 
version  from  other  land  uses  to  Commercial-Industrial  use.  If  the  code  of  an  obser¬ 
vation  has  the  value  of  0.0000,  it  indicates  that  that  pixel  has  been  converted  from 
other  land  use  to  Commercial-Industrial  use.  When  the  code  has  a  value  of  1.0000, 
it  indicates  that  that  pixel  either  had  not  been  changed  or  had  been  changed  to 
Residential  or  Open  Space  use.  According  to  Section  5  of  the  guide  file,  the  order  of 
the  digital  factor  maps  is  “RoadAttCom.asc,”  “SlopeCom.asc,”  and  “WaterAtt- 
Com.asc.”  Therefore,  from  the  second  to  fourth  (1  +  p,  p=3)  columns,  the  values  are 
from  the  factors  represented  by  “RoadAttCom,”  “SlopeCom,”  and  “WaterAttCom,” 
respectively.  Figure  2  illustrates  the  structure  of  the  sampled  data  files. 
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1.0000 

1.0000 

.6269 

.5000 

.0000 

.4499 

.6493 

.5000 

1.0000 

.4499 

.7411 

.5000 

1.0000 

1.0000 

.7513 

.5000 

1.0000 

1.0000 

.5885 

.5000 

Figure  2.  Structure  of  sampled  data  files;  illustrated  file  is  “CI_Sampled.dat”. 
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6  SAS  Program  Builder 


SAS  programs  are  needed  to  estimate  the  weights  (coefficients)  of  LEAM.  The  re¬ 
quired  SAS  programs  are  generated  by  running  “LEAM_SAS_Builder.exe”  (or  first 
compiling  “LEAM_SAS_Builder.f  ’  to  obtain  the  executable  file).  The  generated  SAS 
programs  are  stored  in  subfolder  “Estimatn”. 

Four  SAS  programs  are  generated  by  running  the  SAS  program  builder.  For  each 
land-use  conversion,  one  SAS  program  is  generated  for  estimating  the  weights  of 
the  corresponding  probability  model.  For  the  probability  of  converting  to  land  use 
for  Residential,  Commercial-Industrial,  and  Open  Space,  the  names  of  the  corre¬ 
sponding  SAS  programs  are  respectively  “RES_Est.sas”,  “CI_Est.sas”,  and 
“OS_Est.sas”.  SAS  program  “NCH_Est.sas”  is  for  estimating  the  weights  of  the 
model  to  predict  the  probability  of  “no  conversion”.  In  each  of  these  SAS  programs, 
the  independent  variables  of  the  logistic  models  are  the  factors  and  their  cross 
products  (to  represent  the  factor  interactions).  The  factors  are  named  as  variables, 
and  their  cross  products  are  named  as  variables  (m=[p-l]xp/2).  The  order  of  the 
corresponding  map  file  list  (Sections  5-7  of  the  guide  file)  is  the  order  of  the  inde¬ 
pendent  variables  X’s,  and  is  printed  at  the  beginning  of  the  SAS  program  as  mes¬ 
sage  lines.  Independent  variables  (U’s)  are  automatically  generated  and  defined  in 
the  SAS  programs. 

Since  no  factor  is  defined  for  modeling  “no  conversion,”  the  independent  variables 
X’s  of  this  model  are  created  using  the  independent  variable  X’s  of  the  models  for 
the  three  types  of  conversion.  Independent  variable  U’s  of  the  model  for  “no  conver¬ 
sion”  are  generated  and  defined  in  the  SAS  programs  in  a  manner  similar  to  the 
variables  for  land-use  conversion  described  in  the  preceding  paragraph. 
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7  Integration  of  Estimates 


The  weights  (coefficients)  of  LEAM  are  estimated  by  using  the  SAS  software  pro¬ 
gram.  The  generated  SAS  programs  “RES_Est.sas”,  “CI_Est.sas”,  “OS_Est.sas”,  and 
“NCH_Est.sas”  estimate  the  weights  for  corresponding  land-use  changes  with  spe¬ 
cific  sampled  data  files  (Table  2).  You  will  need  to  manually  generates  data  files 
that  contain  the  estimates  of  the  weights  (coefficients).  Table  2  lists  the  names  of 
those  estimate  files. 

After  running  one  of  the  SAS  programs  listed  in  Table  2,  the  outcome  of  logistic  re¬ 
gression  will  be  displayed  on  the  SAS  Output  Screen  (see  Figure  3  for  an  example). 
Highlight  the  estimates  part  emphasized  in  Figure  3,  copy-paste  that  part  to  Note¬ 
pad,  add  a  Title  Line  and  a  space  line  at  the  top  of  the  context,  then  save  it  as  a  text 
file  with  the  name  according  to  the  SAS  program  and  Table  2.  Figure  4  illustrates 
the  format  of  the  files  of  estimates. 

When  all  files  of  estimates  listed  in  Table  2  have  been  developed,  the  estimates  of 
all  weights  are  integrated  into  one  file,  “Est_All.dat”.  The  executable  file 
“LEAM_Est_Integratn.exe”  (which  can  also  be  obtained  by  compiling 
“LEAM_Est_Integratn.f’)  integrates  the  estimates  of  the  weights  of  LEAM  from  the 
four  independent  files  of  estimates  to  create  “Est_All.dat”.  This  integrated  file  of  es¬ 
timates  is  the  input  source  of  weights  of  LEAM  for  simulation  using  programs  in 
C/C++,  FORTRAN,  or  other  computer  languages.  Figure  5  illustrates  the  format  of 
“Est_All.dat.” 


Table  2.  Types  of  land-use  conversion  and  the  names  of  corresponding  SAS  programs,  sampled 
data  files,  and  estimated  weights  files. 


Land-use  Conversion 

SAS  Program 

Sampled  Data  File 

File  of  Estimates 

Residential  (RES) 

RES_Est.sas 

RES_Sampled.txt 

Est_RES.txt 

Commercial-Industrial  (Cl) 

CI_Est.sas 

CI_Sampled.txt 

Est_Cl.txt 

Open  Space  (OS) 

OS_Est.sas 

OS_Sampled.txt 

Est_OS.txt 

No  Change  (NCH) 

NCH_Est.sas 

Nch_Sampled.txt 

Est_NCH.txt 
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_J  C:\Program  Files\SAS  Institu 


Target 

Block 


Figure  3.  SAS  Output  Screen  displaying  estimates  of  weights  of  LEAM  (for  Residential,  eight  factors 
are  considered).  The  part  labeled  “Target  Block,”  which  is  at  almost  the  end  of  the  Output  Screen, 
is  what  needs  to  be  highlighted,  copied,  and  pasted  into  Notepad  to  create  a  file  of  estimates. 
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Weights  of  LEAM  for  prediction  of  the  probability  of  No  change  <■ 

< - 

Analysis  of  Maximum  Likelihood  Estimates 
Standard  Wald 


Parameter 

DF 

Estimate 

Error 

Chi-Square 

Pr  >  ChiSq 

x2 

1 

1 9 . 5595 

3.9388 

24.6593 

<.0001 

x3 

1 

9 , 8297 

3.4374 

8.1774 

0.0042 

x5 

1 

-34.6420 

4.1868 

68.4593 

<.0001 

x6 

1 

8 , 6733 

1 . 5526 

31  . 2056 

<•0001 

x8 

1 

1  5 . 41  32 

4.1809 

13.5908 

0.0002 

u2 

1 

2.  61  82 

1 .1498 

5.1850 

0.0228 

u3 

1 

3.  7806 

1  . 751  0 

4.6619 

0.0308 

u4 

1 

12.  4261 

2.0478 

36.8205 

<.0001 

u7 

1 

5.  2639 

1 . 8742 

7.8880 

0.0050 

u8 

1 

-3. 4480 

1 .1476 

9.0273 

0.0027 

ull 

1 

-4. 4063 

1 . 5246 

8.3534 

0.0038 

ul  4 

1 

-2.  81  52 

1 . 3525 

4.3329 

0.0374 

u20 

1 

3.  7414 

1 . 7683 

4.4767 

0.0344 

u26 

1 

-15.4770 

3.2647 

22.4745 

<.0001 

Title  line 
Space  line 


Pasted  from 
SAS  Output 
Screen 


Figure  4.  Format  of  the  files  of  estimates.  The  file  of  estimates  in  this  example  is  “Est_NCH.txt,” 
which  is  for  “no  change”  (eight  factors  are  considered). 


Attractor 

Estimates  of 

RES 

the  weights  of  LEAM 

Com_Ind  Open  Space  No_ 

< - 

f - 

Change 

< - 

XI 

- . 368640E+02 

- . 268640E+02  - 

373906E+02 

< - 

. OOOOOOE+OO 

X  2 

. 000000E+00 

- .1  57995E+02 

1 95595E+02 

. 1 95595E+02 

X  3 

. 000000E+00 

,22731  6E+02 

982970E+01 

. 98297 0E+ 01 

X  4 

- . 197995E+02 

. 000000E+00 

OOOOOOE+OO 

, OOOOOOE+OO 

X  5 

.19731  6E  +  02 

. 000000E+00  - 

346420E+02 

- . 34642 0E+ 02 

X  6 

. 000000E+00 

. 607620E+01 

OOOOOOE+OO 

. 86733 0E+ 01 

X7 

. 507620E+01 

. 000000E+00 

867330E+01 

, OOOOOOE+OO 

X  8 

. 000000E+00 

. OOOOOOE+OO 

1 541 32E+02 

.  1  541  32E+02 

Ul 

.10861 5E+02 

.1  1861  5E+02 

OOOOOOE+OO 

. OOOOOOE+OO 

U26 

. 000000E+00 

.000000E+00 

OOOOOOE+OO 

- . 154770E+02 

U27 

. 000000E+00 

- . 234390 E+ 01 

OOOOOOE+OO 

.OOOOOOE+OO 

U28 

. 000000E+00 

•OOOOOOE+OO 

OOOOOOE+OO 

. OOOOOOE+OO 

Title  line 
Space  line 

Space  line 
Space  line 


Figure  5.  Format  of  the  integrated  estimates  (file  “Est_All.dat”).  In  this  example,  eight  factors  are 
considered  in  LEAM. 
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8  Conclusions 


This  report  documents  methodology  and  a  tool  set  to  calibrate  the  driver  weights  for 
LEAM  that  (1)  uses  readily  available  land-use  data,  (2)  requires  minimal  data 
preparation,  (3)  is  easily  implemented,  and  (4)  is  adaptable  for  various  implementa¬ 
tions  of  LEAM.  Implementing  this  methodology  using  the  tools  documented  in  this 
report  provides  optimum  driver  weight  calibrations  based  on  the  available  calibra¬ 
tion  data. 
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Appendix  A:  Estimation  of  the  Weights 

(Coefficients)  of  LEAM 


Based  on  an  earlier  study,  a  logistic  model  including  both  individual  attractors  and 
their  interaction  (product  terms)  was  used  to  predict  the  probability  of  Residential 
development: 

Equation  A1.  Initial  model  for  predicting  the  probability  of  residential  development. 


36 


ae  p(C  I  ton 

t  C  V  k  I  Non-  >  Ci.  /  _i_  o  *  |  o  t  t  *  i  ^  ^  a 

L°g|  - a  A;Xi+a  BjUj  ,  k=l,  2,  3,  4; 

sB"  At'k  L.>cki0  k=  1  j=l 


where  X’s  and  U’s  are  the  attractors  of  LEAM,  their  combinations  listed  in  Table 
Al,  and  A’s  and  B’s  are  weights  (coefficients)  estimated  based  on  the  historical 
maps.  P(  Ck  Res|iVo77  ->Ck)  is  the  probability  of  the  kth  category  of  development. 

When  some  of  the  factors  listed  in  Table  Al  are  missing  in  a  specific  study  area,  the 
corresponding  terms  in  Eq.  Al  will  be  dropped  in  model  calibration  as  well  as  in 
probability  prediction. 

Multinomial  probabilities  were  predicted  in  two  steps.  First,  a  set  of  multinomial 
probability  models  had  been  built  using  selective  logistic  regression  analysis.  Then, 
a  standardization  procedure  was  applied  to  compute  the  multinomial  probabilities 
of  the  categories.  The  following  equations  are  provided  for  each  probability  esti¬ 
mate. 

Equation  A2.  Initial  model  for  predicting  the  probability  of  Commercial-Industrial. 

ae  p'  6  13  36 

A_CI,  XX_CI;  +  a  B  CIj  xU  CIj 

!1-  r  CI 0  i=1  j=1 

Equation  A3.  Initial  model  for  predicting  the  probability  of  Open  Space. 

ae  p'  q  13  3® 

LogS — ^-±=  a  a_0S;  xx_os,  +  a  b  os,  xu_os, 

ii-  p'OS0  i=i  j=i 
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Equation  A4.  Initial  model  for  predicting  the  probability  of  Residential. 

ae  p'  p  13  36 

Logj»  RE  ±=  a  A  RE.  xX  RE  +  a  B  RE.xU  RE. 

Bl  —  id  _  J  —  j 

81 '  rRE0  i=l  j=l 

Equation  A5.  Initial  model  for  predicting  the  probability  of  No  Change. 

ae  p  p  13  36 

Log£  NC  ±=  a  A_NCixX_NCi+  a  B_NCjxU_NCj 
ll-  Pnc0  .-i  j=l 

where  A_Ci  is  the  coefficient  of  the  score  of  attractor  i,  where  C={RE,  Cl,  OS, 
NC}={Residential,  Commercial/industrial,  Open  Space,  No  Change},  X_Wi  is  the 

score  of  attractor  i,  B _ Cj  is  the  coefficient  of  the  product  of  the  scores  of  a  pair  of 

attractors  (listed  in  Table  Al),  U_Cj=Cm*  Cn,  and  NCk=(3-X_REk-X_CIk-X_OSk)/3, 
k=l,  ...,  8. 

To  make  the  predicted  probabilities  consistent,  the  original  probabilities  predicted 
using  the  logistic  models  from  Eqs.  A2-A5  need  to  be  adjusted: 

Equation  A6.  Probability  adjustment  for  initial  model  probabilities. 


P  = 

rLUC 


P' 


LUC 


P'  +  P'  +  P'  +  P' 

1  Cl  1  OS  1  RE  1  NC 


where  LUC={RE,  Cl,  OS,  NC}={land  use  converted  to  Residential,  Commer¬ 
cial/industrial,  Open  Space,  No  Change}. 
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Table  A1.  Codes  of  factors  and  their  combinations  in  Equation  A1. 


Factor 

Code 

Factor 

interaction 

Code 

Factor 

interaction 

Code 

City  Attractor 

XI 

X1*X2 

U1 

X4*X5 

U19 

City  Road  Attractor 

X2 

X1*X3 

U2 

X4*X6 

U20 

Forest  Attractor 

X3 

X1*X4 

U3 

X4*X7 

U21 

Slope  Attractor 

X4 

X1*X5 

U4 

X4*X8 

U22 

Ramp  Attractor 

X5 

X1*X6 

U5 

X5*X6 

U23 

Road  Intersection 

X6 

X1*X7 

U6 

X5*X7 

U24 

State  Highway 

X7 

X1*X8 

U7 

X5*X8 

U25 

Water  Attractor 

X8 

X2*X3 

U8 

X6*X7 

U26 

Neighbor 

X9 

X2*X4 

U9 

X6*X8 

U27 

Utilities 

X10 

X2*X5 

U10 

X7*X8 

U28 

Growth  Trend 

XII 

X2*X6 

U11 

X9*X1 

U29 

Agriculture  Protection 

XI 2 

X2*X7 

U12 

X9*X2 

U30 

Growth  Booster 

X13 

X2*X8 

U13 

X9*X3 

U31 

~~ 

— 

X3*X4 

U14 

X9*X4 

U32 

X3*X5 

U15 

X9*X5 

U33 

X3*X6 

U16 

X9*X6 

U34 

X3*X7 

U17 

X9*X7 

U35 

X3*X8 

U18 

X9*X8 

U36 

Nb_nc  =  0.33*(3-nbR  -  nbC  -  nbO), 

where  {nbR,  nbC,  nbO}  =  {neighbor  Residential,  neighbor  Com.-lnd.,  neighbor  Open_Space}. 
Ut_nc  =  0.33*(3  -  Ut_R  -  Ut_C  -  Ut_0). 
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Appendix  B:  PEWL  Fort  Benning  Example 

In  this  example,  the  land-use  conversion  from  1980  through  1990  inside  and  around 
Columbus  City,  GA,  Phenix  City,  AL,  and  Fort  Benning  are  considered.  The  land- 
use  maps  from  1980  and  1990  and  the  maps  of  the  considered  factors  (attractors) 
are  included  in  subfolder  “Benning_Example”  inside  the  folder  “LEAM_Calb”.  The 
corresponding  guide  file  “Parameters.dat”,  which  lists  all  the  essential  parameters 
and  file  lists  for  weight  estimation,  is  also  included  in  subfolder  “Ben- 
ning_Example”.  In  order  to  run  the  example,  all  map  files  need  to  be  moved  into 
subfolder  “DataFiles”,  and  the  guide  file  “Parameters.dat”  needs  to  be  moved  to 
folder  “LEAM_Calb”.  Then,  follow  the  instruction  in  this  User’s  Guide  to  estimate 
the  weights  of  BEAM. 
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Appendix  C:  FORTRAN  Source  Code 


LEAMSampling.f  Program  Code 

•k’k’k’k’k’k’k’k’k'k’k'k’k’k’k’k’k’k’k’k’k’k'k’k'k’k’k’k’k’k’k’k’k'k’k'k’k’k’k’k’k’k’k’k’k’k'k’k'k’k'k’k’k’k'k’k'k’k'k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k 

c  Programmer:  Shoufan  Fang 

c 

c  Program  to  sample  obs.  for  estimating  the  weights  of  LEAM 
c 

program  Sampler 
c 

USE  MSIMSL 

character  chtmp*100,  FName*50,  filename*40,dirc*85 
integer  nrows,  ncols,  N_S,  N_files,  L_Dirc 

CALL  DISCLMO 

open(unit=30,file="Parameters.dat",  form='formatted') 

read(30,*) 
read(30,*) 
c  Read  Directory 
read(30,*)  dire 

do  i=l,85 

if  (dirc(i:i)  .NE.  ")  L_Dirc=i 
enddo 

read(30,*) 

read(30,*)  chtmp, nrows 
read(30,*)  chtmp, ncols 
c 

read(30,*) 

read(30,*)  chtmp, filename 
FName='DataFiles\'//filename 
chtmp=dirc(l:L_Dirc)//'\LEAM_Calb\'//FName 
open(unit=l,file=chtmp,  form='formatted') 
read(30,*)  chtmp, filename 
FName='DataFiles\'//filename 
chtmp=dirc(l:L_Dirc)//'\LEAM_Calb\'//FName 
open(unit=2,file=chtmp,  form='formatted') 


ERDC/CERL  SR-05-59 


23 


c 


chtmp=dirc(l:L_Dirc)//'\LEAM_Calb\DataFiles\Class_Coords.dat' 

open(unit=15,file=chtmp) 

N_S=0 


CALL  Sampl_LU(nrows,  ncols,N_S) 
CCCC********CCCCCC*******CCCCCCCC 


•k'k'k'k'k'k'k'k 


read(30,*) 
read(30,*)  N_files 


chtmp=dirc(l:L_Dirc)//'\LEAM_Calb\DataFiles\Class_Coords.dat' 
open(unit=25,file=chtmp,  form='formatted') 

CALL  SamplATTR(nrows,ncols,N_S,N_files,dirc) 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  End  of  program 
350  stop 
end 

cccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c  Subroutine  to  determine  which  pixels  should  be  sampled  and  sample  land  use 

conversion 

Subroutine  Sampl_LU(nrs,  ncs,NS) 
c  character  folder*85 

integer  nrs,  ncs,  NS,  n(18) 
real  x(8,nrs,  ncs),a(ncs),rndm 
real*8  acptl,acpt2 

do  nn=l,2 
do  j=l,6 
read(nn,*) 
enddo 
c 

do  i=l,nrs 

read(nn,*)(x(nn,i,j),j=l,ncs) 
c  print*, nn,i,nrows,ncols 
enddo 
close  (nn) 
enddo 

a(l)=1.0*nrs*ncs 
if  (a(l)  .LT.  350000)  then 
acptl=0.025 
else 

acptl=9500.0/(1.0*nrs*ncs) 

endif 
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acpt2=0.5+acptl 

acptl=0.5 


NS=0 

do  i=3,nrs-2 
do  j=3,ncs-2 

a(j)=FUNC_Possbl_LU(x(l,i,j)) 
if  (a(j)  -EQ.  0)  goto  580 
if  (x(l,i,j)  .LE.  -999)  goto  580 

a(j)=  x(l,i,j) 

n(l)=  int(FUN  C_D  vlp_LU  (a(j))) 
a(j)=  x(2,i,j) 

n(2)=  int(FUN  C_D  vlp_LU  (a(j))) 
c 

n(8)=-10 

c 

c  Changed? 

n(3)  =  n(2)-n(l) 
if  (n(3)  .NE.  n(2))  n(3)=0 
c 

rndm=RNUNFO 

if  ((rndm  .GE.  acptl)  .AND.  (rndm  .LE.  acpt2))  then 
NS=NS+1 

ccccccccccc  Initial  land  use 
n(5)=int(x(l,i,j)) 
write(15,201)n(3),n(5),i,j 
endif 

c  read(nn,*)(x(nn,i,j),j=l,ncols) 

580  enddo 
enddo 
close  (15) 

print*, 'Number  of  pixels  sampled:  ',NS 
print*,’ ' 
print*,' ' 

201  format(5(lx,I5)) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c  Sampling  the  factor  maps  according  to  the  determined  coordinates 
Subroutine  SamplATTR(nrs, ncs, NS, N_f, folder) 

character  FName*45,  filename*30,nm_f(3)*15,folder*85,chtmp*100 
integer  nrs,  ncs,  n_cd(4,NS),NS,  NN_S 
real  x(N_f,nrs,ncs) 

data  nm_f  /'RES_sampled.dat', 'CI_sampled.dat', 'OS_sampled.dat'/ 
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do  il=l,NS 

read(25,*)(n_cd(k,il),k=l,4) 
enddo 
close  (25) 

do  il=l,85 

if  (folder(il:il)  .NE.  ")  NN_S=il 
enddo 


do  11=1,3 
NN=II+35 

FName='Sampled\7/nm_f(II) 

chtmp=folder(l:NN_S)//'\LEAM_Calb\7/FName 

open(unit=NN,file=chtmp) 

read(30,*) 

do  LL=l,N_f 

read(30,*)  filename 
FName='DataFilesV//filename 
chtmp=folder(l:NN_S)//'\LEAM_Calb\7/FName 
open(unit=LL,file=chtmp,  form='formatted') 

do  i=l,6 
read(LL,*) 
enddo 

do  i=l,nrs 

read(LL,  *)  (x(LL,i,  j) ,  j = 1 ,  ncs) 
enddo 
close  (LL) 
enddo 


do  i=l,NS 
xy=1.0 

if  (n_cd(l,i)  .EQ.  II)  xy=0.0 

write  (NN,  1 1 0)  xy ,  (x(k,  n_cd(3 ,  i) ,  n_c  d(4  ,i)) ,  k=  1 ,  N_f) 
enddo 
close  (NN) 

enddo 

cccccccccccccccccccccccccccc 

c 

do  11=1,3 

FName='Sampled\7/nm_f(II) 
chtmp=folder(l:NN_S)//'\LEAM_Calb\7/FName 
open(unit=II,file=chtmp,  form='formatted') 
enddo 
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chtmp=folder(l:NN_S)//'\LEAM_Calb\Sampled\Nch_Sampled.dat' 

open(unit=15,file=chtmp) 

do  i=l,NS 
xy=1.0 

if  (n_cd(l,i)  .EQ.  0)  xy=0.0 
do  11=1,3 

read(II,*)xz,(x(k,II,l),k=l,N_f) 

enddo 

Do  LK=1,  N_f 

x(LK,l,5)=(3.0-x(LK,l,l)-x(LK,2,l)-x(LK,3,l))/3.0 

enddo 

write(15,110)xy,(x(k,l,5),k=l,N_f) 

enddo 

do  11=1,3 
close  (II) 
enddo 
close  (15) 


cccccccccccccccccccccccccccc 
110  format(17(lx,f6.4)) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c  Information  subroutine 
Subroutine  DISCLM() 


print*  '*********************************************************' 

print*,'*  *' 

print*,'*  DISCLAIM  *' 

print*,'*  *' 

print*,'*  *' 

print*,'*  The  authors  of  this  program  are  not  responsible  *' 

print*,'*  to  any  computer  hardware  damage  or  information  loss  *' 

print*,'*  due  to  improper  use  of  this  program.  *' 

print*,'*  *' 

print*,'*  *' 

print*,'*  *' 

print*,'*  ALL  RIGHTS  OL  THIS  PROGRAM  ARE  RESERVED 
print*,'*  *' 

print*,'*  *' 

print*  I********************************************************** 
print*,'  ' 
print*,'  ' 


ALL  RIGHTS  OF  THIS  PROGRAM  ARE  RESERVED 


pause 
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c 


c 


do  LK=1,10 
print*,' ' 
enddo 


print*  '*********************************************************' 
print*,' 

print*,'  ' 

print*,'  Sampling  for  estimation  of  the  weights  of  LEAM  ' 

print*,' 

print*,'  ' 

print*  '*********************************************************' 


print*,'  ' 
do  LK=1,5 


print*,' ' 
enddo 


return 

end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
function  FUNC_WATER(w) 

real  w,  q 


q  =  0 

IF  ((w  .EQ.  11)  .or.  (w  .EQ.  12))  q  =  1 

FUN  C_WATER=q 

return 

end 

ccc******+++++ . ===== 

c  The  folowings  are  working  functions  according  to  LEAM 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
function  FUNC_CMI(w) 

real  w,  q 


q  =  0 

IF  ((w  .GE.  23)  .and.  (w  .LE.  24))  q  =  1 

FUNC_CMI=q 

return 

end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
function  FUNC_OPN(w) 

real  w,  q 


q  =  0 

IF  (w  .EQ.  85)  q  =  1 

FUNC_OPN=q 

return 
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end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
function  FUNC_RES(w) 

real  w,  q 

q  =  0 

IF  ((w  .GE.  21)  .and.  (w  .LE.  22))  q  =  1 

FUNC_RES=q 

return 

end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
function  FUNC_Dvlp_LU(w) 

real  w,  q 

q  =  0 

IF  ((w  .EQ.  21)  .or.  (w  .EQ.  22))  q  =  1 
IF  ((w  .EQ.  23)  .or.  (w  .EQ.  24))  q  =  2 
IF  (w  .EQ.  85)  q  =  3 

FUN  C_Dvlp_LU=q 

return 

end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
function  FUNC_Possbl_LU(w) 

real  w,  q 

q  =  1 

IF  ((w  .EQ.  21)  .or.  (w  .EQ.  22))  q  =  0 
IF  ((w  .EQ.  23)  .or.  (w  .EQ.  24))  q  =  0 
IF  ((w  .EQ.  11)  .or.  (w  .EQ.  85))  q  =  0 

FUNC_Possbl_LU=q 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccc 


•k'k'k'k'k'k'k'k'k'k’k’k'k'k'k'k'k'k'k'k'k'k'k’k’k'k'k'k'k'k'k'k'k'k’k'k’k'k'k'k'k'k'k'k'k’k’k'k'k’k’k'k'k’k'k'k'k’k’k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 
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LEAM  SAS  Builder.f  Program  Code 

'k'k'k'k'k'k'k’k’k’k’k’k’k'k'k'k'k'k'k'k’k’k’k’k’k'k'k'k'k'k'k’k’k’k’k’k’k'k'k’k'k'k'k’k’k’k’k'k’k’k’k'k’k’k’k'k’k’k’k'k’k’k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k’k 

c  Programmer:  Shoufan  Fang 

c 

c  Program  to  generate  SAS  programs  for  estimating  the  weights  of  LEAM  model 
c 
c 
c 

program  SAS_Builder 
c 
c 

Character  folder*85 
Integer  N_files 

CALL  DISCLMO 

open(unit=30,file="Parameters.dat",  form='formatted') 


Do  i=l,  10 

if  (i  .EQ.  3)  then 
read(30,*)  folder 
else 

read(30,*) 

endif 

Enddo 

read(30,*)  N_files 
read(30,*) 
c 

CALL  SAS_PROGRM(N_files, folder) 
c 

close  (30) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  End  of  program 
350  stop 
end 

cccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
Subroutine  SASJNTI(M) 

Integer  M 

Write  (M,*)  'options  ps=1000  ls=90;' 

Write  (M ,*) 


return 


30 


ERDC/CERL  SR-05-59 


end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
Subroutine  SAS_PROGRM(N_F,fldr) 
integer  N_F,  N_VAR,  nn(10),  N_FL 

character  chtmp*125,  VNa m e (N_F) * 5 0 ,  SAS_NM(4)*12,cr(2)*l 
character  F_nm(4)*15,DtNm(4)*5,resp(4)*3,cc(5)*3  ,fldr*85 
c 

data  F_nm/'RES_Sampled. datVCI_Sampled.datVOS_Sampled.dat', 

+  'NCh_Sampled.dat7 
data  DtNm/'one','two','three','four7 
data  resp/'RES','Cr,'OS','NCh7 

data  SAS_NM/'RES_Est.sas','CI_Est.sas','0S_Est.sas','NCh_Est.sas7 
c 

do  i=l,85 

if  (fldr(i:i)  .NE.  ")  N_FL=i 
enddo 

do  i=l,4 

chtmp=fldr(l:N_FL)//"\LEAM_Calb\Estimatn\'7/SAS_NM(i) 

open(unit=i,file=chtmp) 

CALL  SAS_INTI(i) 
enddo 

cr(l)=CHAR(39) 

cr(2)=CHAR(59) 

do  LL=1,N_F 

CALL  N_CHAR(LL,  cc(l)) 
read(30,*)  VName(LL) 
do  JK=l,len(VName(LL)) 

if  (VName(LL)(JK:JK)  .NE.  ")  nn(10)=JK 
enddo 

nn(10)=nn(10)+l 

chtmp="*  "//"  x'7/cc(l)//"=  '7/VName(LL)(l:nn(10))//cr(2) 
do  i=l,4 

Write  (i,*)  chtmp 
enddo 
enddo 
do  i=l,4 

Write  (i,*) 

Write  (i,*) 
enddo 

c  F_nm-_Sample.dat' 

N_VAR=N_F*(N_F- 1)/2 

CALL  N_CHAR(N_F,  cc(5)) 


CALL  N_CHAR(N_VAR,  cc(4)) 
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do  JJ=1,4 

chtmp="Data  "//DtNm(JJ)//cr(2) 

Write  (JJ,*)  chtmp 

CALL  GETFILE(F_nm(JJ),resp(JJ),cc(5),JJ,fldr) 
nn(JJ)=0 
do  11=1,  N_F  -1 
do  IJ  =  II+l,  N_F 
nn(JJ)=nn(JJ)+l 
CALL  N_CHAR(nn(JJ),  cc(l)) 

CALL  N_CHAR(II,  cc(2)) 

CALL  N_CHAR(IJ,  cc(3)) 
chtmp="U"//cc(l)//"=  x"//cc(2 )//"  *  x"//cc(3)//cr(2) 

Write  (JJ,*)  chtmp 
enddo 
enddo 

Write  (JJ,*) 

Write  (JJ,*) 

CALL  SETMODEL(resp(JJ),DtNm(JJ),cc(5),cc(4),JJ) 
close  (JJ) 
enddo 


c  201  format(5(lx,I5)) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  GETFILE(Fname,rsp,cnv,M,fold) 
character  Fname*15,chr*100,cnv*3,cr(2)*3,rsp*3,fold*85 
Integer  M  ,  NFD 

do  KL=1,85 

if  (fold(KL:KL)  .NE.  ")  NFD=KL 
enddo 

cr(l)=CHAR(39) 

cr(2)=CHAR(59) 

chr="infile  "//cr(l)//fold(l:NFD)//"\LEAM_Calb\Sampled\"// 

+  Fname//cr(l)//cr(2) 

Write  (M,*)  chr 

chr="Input  "//rsp//"  xl  -  x"//cnv//cr(2) 

Write  (M,*)  chr 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  N_CHAR(n5,  ccc) 
character  ccc*3,  ct(2)*l 
Integer  n5,  nl,  n2 
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nl=int(n5/10) 

n2=n5-10*nl 

ct(l)=CHAR(48+n2) 

ccc=ct(l) 

if  (nl  .GT.  0)  then 
ct(2)=CHAR(48+nl) 
ccc=ct(2)//ct(l) 
endif 


return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  SETMODEL(rsp,DataNm,cvn,ccvn,M) 
character  chr*85,DataNm*5,rsp*3,ccr*3,cvn*3,ccvn*3 
Integer  M 

ccr=CHAR(59) 

chr="PROC  LOGISTIC  data=7/DataNm//ccr 
Write  (M,*)  chr 

chr="Model  "//rsp//"  =  xl  -  x"//cvn//"  U1  -  U"//ccvn// 

+  "  /  noint  rsq  selection=backward  "//ccr 

Write  (M,*)  chr 

chr="RUN"//ccr 

Write  (M,*)  chr 

Write  (M,*) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  DISCLM() 


print*, 

4-* 


• k  I'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'ki 

*1 

DISCLAIM  *' 


print*, 
print*,'* 

print*,'*  *' 

print*,'*  *' 

print*,'*  The  authors  of  this  program  are  not  responsible  *' 
print*,'*  to  any  computer  hardware  damage  or  information  loss  *' 
print*,'*  due  to  improper  use  of  this  program.  *' 

print*,'* 
print*,'* 
print*,'* 
print*,'* 
print*,'* 
print*,'* 

print*  '*********************************************************' 
print*,’  ’ 
print*,’  ’ 


** 

** 

ALL  RIGHTS  OF  THIS  PROGRAM  ARE  RESERVED  *’ 

** 


pause 
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c 


do  LK=1,10 
print*,' ' 
enddo 


print*  '*********************************************************' 
print*,' 
print*,' 
print*,' 
print*,' 
print*,' 

print*,'  ' 

print*  '*********************************************************' 


Generating  SAS  Program  for 

f 

estimation  of  the  weights  of  LEAM 


print*,' ' 
do  LK=1,5 


print*,' ' 
enddo 


return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

•k’k’k’k’k’k’k’k’k’k’k'k’k’k’k’k’k’k’k’k’k’k’k’k'k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k'k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k’k 


LEAM  Est  lntegratn.f  Program  Code 

'k'k’k'k’k'k'k’k’k’k’k’k’k'k'k'k'k’k'k'k’k’k’k’k’k'k'k'k'k'k'k’k’k’k’k’k’k'k'k'k'k'k'k’k’k'k’k’k’k’k’k’k’k'k’k'k’k’k’k'k’k’k'k'k'k’k'k’k'k'k'k’k'k’k'k'k'k’k 

c  Programmer:  Shoufan  Fang 

c 

c  Program  to  integrate  all  estimates  of  the  weights  of  LEAM 
c 
c 
c 

program  Est_integrator 
c 
c 

Integer  N_files  ,  N_V 
CALL  DISCLMO 

open(unit=30,file="Parameters.dat",  form='formatted') 

Do  i=l,  10 
read(30,*) 

Enddo 


read(30,*)  N_files 
read(30,*) 
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c 

close  (30) 
c 

N_V=N_files*(N_files- 1)/2 
print*,  N_files,  N_V 
CALL  Integrtn(N_files,  N_V) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  End  of  program 
350  stop 
end 

cccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
Subroutine  OutFile() 
character  chtmp*125 


open(unit=15,file='Estimatn\Est_All.dat') 

Write  (15,*) '  Estimates  of  the  weights  of  LEAM' 

Write  (15,*) 

chtmp-  Attractor  RES  Com_Ind  Open  Space  No_Change' 
Write  (15,*)  chtmp 
Write  (15,*) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
Subroutine  Integrtn(N_F,  N_VAR) 
integer  N_F,  N_VAR,  nn(10) 

character  chtmp*55,  SAS_NM(4)*12,crr*4,vcode*5 
real*8  x(4,N_F),  U(4,N_VAR),a(10) 
c 

data  SAS_NM/'Est_RES.txt','Est_CLtxt7Est_OS.txt7Est_NCh.txt7 
c 

CALL  OutFileQ 


do  i=l,4 

do  j=l,N_F 
x(i,j)=0 
enddo 

do  j=l,N_VAR 


U(i,j)=0 


enddo 

enddo 
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do  235  i=l,4 

chtmp="  Estimatn\"//SAS_NM(i) 
open(unit=i,file=chtmp,  form='formatted') 
do  LL=1,7 
Read  (i,*) 
enddo 

112  read(i,*,  end=225)  err,  (a(J),J=l,3) 

CALL  Decode(crr,nn(l),nn(2)) 
if  (nn(2)  .EQ.  1)  then 
x(i,nn(l))=a(2) 
else 

U(i,nn(l))=a(2) 
endif 
goto  112 

225  close(i) 

235  continue 


do  LL=1,N_L 

CALL  N_CHAR(LL,  err) 
vcode='X7/crr 

Write  (15,350)  vcode,(x(j,LL),j=l,4) 
enddo 

do  LL=1,N_VAR 

CALL  N_CHAR(LL,  err) 
vcode='U7/crr 

Write  (15,350)  vcode,(U(j,LL),j=l,4) 
enddo 

350  format(2x,A5,3x,5(E12.6,2x)) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  Decode(ccc,nl,n2) 
character  ccc*4,  ct(4)*l 
Integer  n5,  nl,  n2,  n3,  nm(5) 

n5=len(ccc) 


do  LL=l,n5 

ct(LL)=ccc(LL:LL) 
if  (ct(LL)  .NE.  ")  n3=LL 
enddo 

nl=0 
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do  LL=2,n3 

nm(5)=10**(n3-LL) 

nm(LL)=ichar(ct(LL))-48 

nl=nl+nm(LL)*nm(5) 

enddo 


n2=2 

if  (ct(l)  .EQ.  'x')  n2=l 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  N_CHAR(n5,  ccc) 

character  ccc*4,  ct(2)*l 

Integer  n5,  nl,  n2 

nl=int(n5/10) 

n2=n5-10*nl 

ct(l)=CHAR(48+n2) 

ccc=ct(l) 

if  (nl  .GT.  0)  then 
ct(2)=CHAR(48+nl) 
ccc=ct(2)//ct(l) 
endif 


return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  DISCLM() 


print*, 

4-* 


• k  I'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'ki 

*1 

DISCLAIM  *' 


print*, 
print*,'* 

print*,'*  *' 

print*,'*  *' 

print*,'*  The  authors  of  this  program  are  not  responsible  *' 
print*,'*  to  any  computer  hardware  damage  or  information  loss  *' 
print*,'*  due  to  improper  use  of  this  program.  *' 

print*,'* 
print*,'* 
print*,'* 
print*,'* 
print*,'* 
print*,'* 

print*  '*********************************************************' 
print*,’  ’ 
print*,’  ’ 


** 

** 

ALL  RIGHTS  OF  THIS  PROGRAM  ARE  RESERVED  *’ 

** 


pause 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


do  LK=1,10 
print*,' ' 
enddo 


print*  '*********************************************************' 
print*,' 
print*,' 
print*,' 
print*,' 
print*,' 

print*,'  ' 

print*  '*********************************************************' 
print*,'  ' 
do  LK=1,5 
print*,' ' 
enddo 


Integrating  the  estimation 

! 

of  the  weights  of  LEAM 


return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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